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Scalar-tensor theories are the best motivated alternatives to general relativity and provide a 
mathematically consistent framework to test the various observable predictions. They can 
involve three functions of the scalar field: (i) a potential (as in “quintessence” models), (ii) a 
matter-scalar coupling function (as in “extended quintessence”, where it may also be rewritten 
as a nonminimal coupling of the scalar field to the scalar curvature), and (iii) a coupling 
function of the scalar field to the Gauss-Bonnet topological invariant. We recall the main 
experimental constraints on this class of theories, and underline that solar-system, binary- 
pulsar, and cosmological observations give qualitatively different tests. We finally show that 
the combination of these data is necessary to constrain the existence of a scalar-Gauss-Bonnet 
coupling. 


1 Introduction 

In the most natural alternative theories to general relativity (GR), gravity is mediated not 
only by a (spin-2) graviton corresponding to a metric g fW , but also by a (spin-0) scalar field 
ip. Such scalar partners generically arise in all extra-dimensional theories, and notably in string 
theory. A dilaton is indeed already present in the supermultiplet of the 10-dimensional graviton, 
and several other scalar fields (called the moduli) also appear when performing a Kaluza-Klein 
dimensional reduction to our usual spacetime. They correspond to the components of the metric 
tensor g mn in which m and n label extra dimensions. Moreover, contrary to other alternative 
theories of gravity, scalar-tensor theories respect most of GR’s symmetries: conservation laws, 
constancy of non-gravitational constants, and local Lorentz invariance even if a subsystem is 
influenced by external masses. They can also satisfy exactly the weak equivalence principle 
(universality of free fall of laboratory-size objects) even for a massless scalar field. 
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Scalar fields are also involved in the cosmological models which reproduce most consistently 
present observational data. In particular, inflation theory is based on the presence of a scalar p 
in a potential V{<p) (for instance parabolic). It behaves as a fluid with a positive energy density 
STrGpp = p 2 + 2V(tp) but a negative pressure 8 irGpp = <p 2 — 2V(p). This causes a period of 
exponential expansion of the universe, which can explain why causally disconnected regions at 
present may have been connected long ago. The isotropy of the observed cosmic microwave 
background (CMB) can thus be understood. Inflation also predicts that our universe is almost 
spatially flat, just because any initial curvature has been exponentially reduced by the expansion. 
This is in remarkable agreement with the location of the first acoustic pe ak o f the CMB spectrum 
at a multipolar index^ £ ~ 220. Observations of type la supernovaetell us that there is 
about 70% of negative-pressure dark energy in our present universe (Ha — 0.7), suggesting that 
its expansion has been re-accelerating recently (since redshifts z ~ 1). This can be explained 
by the presence of a cosmological constant A in GR, but the quantity Ha ~ 0.7 translated 
in natural units gives an extremely small value A ~ 3 x 10 ~ 122 c 3 /(HG), very problematic for 
particle physics if A is to be interpreted as the vacuum energy. This is the main reason why 
“quintessence” models have been proposed, in which the cosmological constant is replaced again 
by the potential V(ip) of a scalar field. Its evolution towards a minimum of V during the 
cosmological expansion then explains more naturally why the present value V(po) ~ A/2 is so 
small. 

Besides these theoretical and experimental reasons for studying scalar-tensor theories of 
gravity, one of their greatest interests is to embed GR within a class of mathematically consistent 
alternatives, in order to understand better which theoretical features have been experimentally 
tested, and which can be tested further. 

The following action defines the most general theory satisfying the weak equivalence principle 
and involving only one spin-0 degree of freedom besides the usual (spin-2) graviton: 

S = Matter [matter;^ = A 2 {p)g tiv \ 

-hj V=9 W(p) (R 2 ^ - + R 2 ) , (1) 

where g^ u is the metric to which matter is universally coupled, and g^y is the Einstein metric 
(describing the spin-2 degree of freedom). This action involves three function of the scalar 
field: a coupling function A(ip) to matter, a potential V(p), and a coupling function W(p) to 
the Gauss-Bonnet topological invariant. Any other combination of the curvature tensor would 
introduce an extra scalar field in the theory, and/or a second negative-energy massive graviton 
which would make the model unstable. 

In Sections 2 and 3, we will not consider any scalar-Gauss-Bonnet coupling, and set W(p) = 
0. The potential V(p) will also be neglected in Section 2, in which we will review the main 
experiment al co nstraints on the coupling function A{p), coming from solar-system and binary- 
pulsar dataP^ In Section 3, we will summariz e ou r results concerning the reconstruction of 
A(p) and V(p) from cosmological observations^^ Section 4 will be devoted to the scalar- 
Gauss-Bonnet coupling W(p), which can be constrained only if one takes into account both 
solar-system and cosmological dataP^We will finally give our conclusions in Section 5. 

2 Solar-system and binary-pulsar constraints 

The effects of a massive scalar field has a negligible effect on the motion of celestial bodies if its 
mass is large with respect to the inverse of the interbody distances. On the other hand, if its 
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Figure 1: Solar-system and binary-pulsar constraints on the matter-scalar coupling function ln^4(<p) = 0:0(93 — 
ipo) + |/3o(93 — 930) 2 + 0(ip — ipo) 3 - The allowed region is shaded. The vertical axis (/3o = 0) corresponds to 
Brans-Dicke theory with a parameter 2 wbd + 3 = 1 /o:q . The horizontal axis (oo = 0) corresponds to theories 
which are perturbatively equivalent to GR, i.e., which predict strictly no deviation from it (at any order 1/c”) in 
the weak-field conditions of the solar system. 


mass is small enough, its potential V(<p) can be locally neglected, but its coupling function to 
matter, A(<p), is strongly constrained by experiment. 

The predictions of metric theories of gravity in weak-field conditions can be parametrized by 
a set of 10 real numbers in the so called “PPN” formalism (parametrized post-Newtonian). All 
of them are presently constrained to be very close to their general relativistic values, and in par¬ 
ticular the t wo famous Eddington parameters f3 and 7 (both equal to 1 in GR). In scalar-tensor 
theoriesW^ they are related to the first two derivatives of In A(ip). computed at the background 
value po of the scalar field. They give the constraints displayed as a thin line in Figure [T] (where 
the Moon symbol refers to Lunar Laser Ranging, the Mercury symbol to the perihelion shift of 
this planet, and the star symbol to light deflection as measured by Very Long Baseline Interfer¬ 
ometry). Solar-system tests thus constrain the first derivative ao = 9 In A((p)/dip to be small, 
but do not tell us much about the second derivative [3q = d 2 In A(ip)/dip 2 . If ao is small enough, 
arbitrary large positive or negative values of /?o are a priori allowed. 

Binary-pulsar give qualitatively different constraints because of nonperturbative strong-field 
effects. Indeed, the largest deviations from the flat metric, at the surface of a star, are of or der 
GM/Rc 2 ~ 0.2 for a pulsar (neutron star), as compared to 2 x 10 -6 for the Sun. We showed^El 
that if (3o is negative, then it is energetically favorable for a neutron star, above a critical mass, 
to create a nonvanishing scalar field. Since this is analogous to the spontaneous magnetization of 
ferromagnets, we called this effect “spontaneous scalarization”. Such macroscopic scalar charges 
change drastically the physics of a binary system, notably because it emits dipolar gravitational 
(scalar) waves oc 1 /c 3 , much larger that the usual quadrupolar radiation oc 1 /c 5 predicted by 
GR. This is the reason why binary-pulsar data, which are consistent with GR, rule out scalar- 
tensor models such that /3o < —5, even for vanishingly small values of ao (he., even if they are 
strictly indistinguishable from GR in the solar system). 

We also showed that the LIGO/VIRGO interferometers will be more sensitive to /?o than 
solar-system tests, but binary-pulsar data are so precise that they already exclude the models 
which predict significant effects in the gravitational waveforms. This is a good news, since it 
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proves that pure GR wave templates suffice to analyze future LIGO/VIRGO data. On the other 
hand, it was shown^ that the LISA interferometer can be sensitive to scalar effects which are 
still allowed by all present tests. 

In conclusion, solar-system tests tightly constrain the first derivative of In A( 93 ) (linear 
matter-scalar coupling strength), whereas binary-pulsar data impose that its second deriva¬ 
tive (quadratic coupling matter-scalar-scalar) is not large and negative. We will now see that 
cosmological observations give access to the full shape of this coupling function, of course not 
with the same accuracy as the above tests, but with the capability of constraining any higher 
derivative of In A((p) (vertex of matter with any number of scalar lines). Moreover, cosmological 
data can also give access to the full shape of the potential V(ip). 

3 Reconstruction of a scalar-tensor theory from cosmological observations 

In cosmology, the usual approach to study quintessence models is to assume a particular form for 
the potential V(<p) (and the matter-scalar coupling function A{<p) when one considers “extended 
quintessence” models), to compute all possible observable predictions, and to compare them to 
experimental data. 

In contrast, in the phenomenological approach, one wishes to reconstruct the Lagrangian of 
the theory from cosmological observations. We proved^ that the knowledge of the luminosity 
distance Dl(z) and of the density fluctuations S m (z) = Sp/p as functions of the redshift 2 indeed 
suffices to reconstruct both the potential V{<p) and the coupling function A{<p). Although the 
explicit reconstruction needs some algebra, this result seems anyway obvious: It is possible to 
fit two observed functions [Dl(z) and S m (z)] thanks to two unknown ones [V (<p) and A(p)\. 

However, future experiments (like the SNAP satellite) will only give access to the luminosity 
distance Dl(z ) with a good accuracy, and the density contrast S m (z) cannot yet be used to 
constrain the models. A semi-phenomenological approach can thus be useful: We make some 
theoretical hypotheses on either the potential V(<p) or the coupling function A(tp), and we 
reconstruct the other one from D[,(z). A priori, one may think that such a reconstruction is 
again obvious: We fit one observed function [Dl(z)] with one unknown function [V(ip) or A(<p)\. 
But this naive reasoning is only valid locally, on a small interval. Indeed, the reconstructed 
function may for instance diverge for some value of the redshift, or one of the degrees of freedom 
may need to take a negative energy beyond a given redshift, which would make the theory 
unstable (and ill defined as a field theory on the surface where the energy changes its sign). 
The positivity of the graviton energy implies A 2 (ip) > 0, which can be translated in terms of 
the standard Brans-Dicke scalar field as 3 >bd > 0. On the other hand, the positivity of the 
scalar-field energy imposes the minus sign in front of the scalar kinetic term —(d^ip) 2 in action 
0 , which translates as cubd > — § in terms of the standard Brans-Dicke parameter. We showed 
that these conditions impose tight constraints on the theories as soon as one knows Dl(z ) over 
a wide enough interval z £ [ 0 , ~ 2 ]. 

For instance, we proved that the present accelerated expansion of the universe can be per¬ 
fectly described by a scalar-tensor theory with a vanishing potential V(tp) = 0 (and therefore 
a vanishing cosmological constant too). We derived analytically the coupling function A(<p) 
which reproduces exactly the same evolution of the scale factor a(z) as the one predicted by 
GR plus a cosmological constant. We even found that the reconstructed function In A(<p) has a 
nice parabolic shape, with a minimum very close to the present value <po of the scalar field, and 
a positive second derivative. This is not only consistent with binary-pulsar data (which forbid 
large and negative values of this second derivative) but also with the cosmological attractor phe¬ 
nomenon analyzed by Damour and Nordtvedt^ The scalar field is generically attracted towards 
a minimum of In A(tp) during the cosmological expansion, whereas some fine tuning would be 
necessary to reach a maximum (negative second derivative). Therefore, we are in the difficult 
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situation in which two very different theories are both consistent with experimental data, and 
there seems to be no way to distinguish them. Fortunately, we found that this scalar-tensor 
theory cannot mimic GR plus a cosmological constant beyond a redshift z ~ 0.7, because the 
scalar field ip would diverge at this value, and above all because the graviton energy would be¬ 
come negative beyond. Therefore, it suffices to measure Dl(z) precisely enough up to 2 ~ 1 to 
rule out such a potential-free scalar-tensor theory. Actually, if Dl(z) is measured over a wider 
interval z & [0, ~ 2], we showed that large experimental errors (tens of percent) were not prob¬ 
lematic: It is still possible to distinguish this potential-free model from GR plus a cosmological 
constant, and thereby to rule out one of them. The results of the SNAP satellite up to z ~ 2 
will therefore be very useful to constrain scalar-tensor theories of gravity. 

One can also impose a particular form of the coupling function A{(p) and reconstruct the 
potential V{<p) which reproduces the observed luminosity distance Dl(z). For instance, for a 
minimally coupled scalar field A{tp) = 1 (usual “quintessence”) in a spatially curved universe, 
we analytically derived the expression of V(ip) which gives the same cosmological evolution 
as GR plus a cosmological constant in a spatially flat universe. We found that the shape 
of the potential is smoother when the universe is (marginally) closed. If it is flat or almost 
flat, one obviously recovers a cosmological constant with its unnaturally small value A ~ 3 X 
l()~ 122 c 3 / (hG). Therefore, in that case, aesthetic reasons may help us discriminate between the 
theories, instead of the much stronger argument of the positivity of energy that we used above. 
This shows anyway that the sole knowledge of Dl(z) suffices to constrain scalar-tensor theories 
of gravity. 

4 Experimental constraints on a scalar—Gauss-Bonnet coupling 

In order to illustrate the different kinds of experimental constraints that can be irnposed^^ on 
the scalar-Gauss-Bonnet coupling function W(<p), we will now focus of a theory with A(ip) = 1 
and V ((f) = 0 in action JTJ. 

Solar-system and binary-pulsar tests are local, and any deviation from GR depends on the 
magnitude of the scalar field created by a massive body. Let us thus analyze first the equation 
satisfied by ip in the vicinity of a spherical mass Mq . We can assume that the metric is close to 
the Schwarzschild solution, and we get at the first nonvanishing order in powers of GM 0 /c 2 



where we have set 7 q = \imGhj (? , and where the derivative W\<p) has been expanded in powers 
of <p in the right-hand side. Since we are assuming that <p takes small values, let us neglect the 
contribution W{(<p. We can then compute any observable prediction, but we quote below only 
the results for the light deflection angle (A0*) and for the perihelion shift per orbit (A 9 P ), which 
suffice for our purpose: 


A 0* 

4 GMq 1536 / 

p 0 c 2 + 35 \ 

WP 


(3) 

a e p 

67rGM 0 

= -^ + 192tt 

pc z 

(GMq \ 2 
V pc 2 ) 

(?)*■ 

(4) 


where po is the minimal distance between the light ray and the Sun, and p is the semilatus 
rectum of an orbit. The first terms on the right-hand sides are the usual general relativistic 
predictions, at first order in GM 0 /c 2 . In conclusion, solar-system (and binary-pulsar) tests can 
easily be passed if |Wq| is small enough. 
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Figure 2: Scalar-Gauss-Bonnet coupling function W(ip) which exactly reproduces the cosmological expansion 
predicted by GR plus a cosmological constant. 


One can now reconstruct the full shape of W ( 95 ) from the cosmological observation of the 
luminosity distance Dl(z), as in the previous section. We found that this can always been 
done, without any problem of negative energy, contrary to what we saw in Section 3. Moreover, 
there exists again an attraction mechanism which drives the scalar field towards a minimum of 
W(tp) during the cosmological expansion. Therefore, a small value of the slope |Wq| is indeed 
expected at present, consistently with what is needed for solar-system tests. In conclusion, we 
are faced again with a serious problem: We just found a theory which seems to be consistent 
with all experimental data, although it is very different from GR in its field content. Our aim 
is therefore to find a way to distinguish it from GR, or to rule it out for internal consistency 
reasons. 

Figure [ 2 ] displays this reconstructed coupling function W (ip ), in which the present value of 
the scalar field is close to the minimum at <p = 0. Its shape is nicely smooth, but its second 
derivative at the origin is huge if one divides it by the tiny factor (16irGh/c 5 )HQ , where H 0 
denotes the Hubble constant. One gets Wq ~ 7x 10 119 , which is in fact not surprizing, since the 
coupling function W(<p) behaves in action (JT|) as the inverse of a cosmological constant. Indeed, 
W(ip) multiplies the square of the curvature tensor, whereas the usual Einstein-Hilbert term 
involves the first power of the scalar curvature, and a cosmological constant does not multiply 
any curvature term at all. Therefore, it was expected that W(cp) involve a dimensionless number 
of the order of the inverse of ( hG/c 3 )A ~ 3 x 10 -122 . Therefore, this model is ugly, but it is not 
yet ruled out. One should not confuse fine tuning and large (or small) dimensionless numbers 
in a model. We are here in the second situation, but there is a priori no fine tuning since 
the scalar field is attracted towards the minimum of W (<p) during the cosmological expansion. 
There remains to study how efficiently it is attracted, but this is actually not necessary for our 
purpose. 

Indeed, Wq takes such a gigantic value that an approximation that we made to analyze 
solar-system tests is no longer valid. Indeed, we have |WgVl 2 > | Wq |, so that the second term on 
the right-hand side of Eq. (J2J cannot be neglected. To simplify the discussion, we will anyway 
assume that W(<p) is parabolic, which is a good approximation in a vicinity of the minimum 
ip = 0 . We will thus neglect the higher order terms 0(ip 2 ) in Eq. (J2J) ; taking them into account 
would not change our conclusions below. We did not find a close analytic solution to Eq. 






but is is possible to write it as a series 




M V_i_ 

Wg n>i (3 x 4)(7 x 8) • • ■ (4n — l)(4n) 


Wg 


[ cos 
cosh 


( GMqTq 
\ r 2 c 2 




if Wg < 0, 
if Wg > 0. 


(5) 

( 6 ) 


The second expression is a good approximation if the argument of the cosine (or hyperbolic 
cosine if Wg > 0) is much greater than 1. This is the case if we use the huge value of Wg 
obtained above from the cosmological reconstruction, and a typical solar-system distance for 
the radius r: The argument of the hyperbolic cosine is then of order 10 8 . 

The above solution is such that ip cx Wg therefore we do not find any nonperturbative 
effect similar to the “spontaneous scalarization” of neutron stars mentioned in Section 2 above. 
Moreover, ^->0asr->oo, and we recover GR for distances r > 4 x 10 14 m (i.e., farther that 
the solar system including Oort’s comet cloud). On the other hand, there are highly nonlinear 
corrections proportional to l/r 4n within the solar system. Since the ratio (12r ( j G 2 M| Wg/r 4 c A ) 
is much greater that 1, its successive powers blow up, but they are compensated by the factors 
1/(3 x 4 x 7 x • 4n) which behave like the inverse of factorials. Therefore, the successive terms 
of series © start to grow exponentially, then reach a maximum for a value of the index n which 
may be large, and finally tend towards zero. Each of these successive terms must be assumed 
to be small enough for the model to pass all classical tests, but one should not forget that the 
largest one does not correspond to n = 1. 

In order to study the effects of such highly nonlinear terms in the solar system, we compute 
their corrections to the Schwarzschild metric in the form 


ds2 - - 1 +^ ° 2dt2 + 1 + 5Z \ d p 2 + 


f- sin 2 6 d(j) 2 ) 


and we find that the light deflection angle and the perihelion shift are respectively given by 


) ^n(n-l)/3 n c 2 2 

2OA/0P- 1 " + 


(7) 

Y 

( 8 ) 
(9) 


which generalize Eqs. ©-© above. Note that these results are (independently) perturbative 
in each of the coefficients a n and (3 n , because we are assuming that the scalar-field effects are 
negligible with respect to the general relativistic predictions. However, the dominant scalar-field 
corrections may correspond to a large value of index n. 

When solution © or its approximation © are used to compute the metric coefficients a n 
and p n in the above observable predictions, and if we use the huge value of Wg obtained from 
the previous cosmological reconstruction of W(ip), we get the following experimental constraint: 

|Wo| < 10 _2xl ° n . (10) 


Now we can speak of fine tuning, and even of hyperfine tuning! This constraint simply means 
that the present value of the scalar field must be exactly at the minimum of the coupling 
function W{tp), otherwise solar-system tests are violated. And since the universe is still evolving, 
the scalar field cannot remain so close to the minimum for more than a fraction of a second. 
Therefore, even if we assumed that W/ = 0 strictly to pass solar-system tests, this would not 
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be the case a tiny instant later. In conclusion, we managed to rule out the scalar-tensor model 
A{<p) = 1 , V(<p) = 0 and W{p) ^ 0. It cannot describe an accelerating expansion of the universe 
at present and pass solar-system (and binary-pulsar) tests at the same time. 

Of course, this result does not rule out any scalar-Gauss-Bonnet coupling. A model with 
three (or even two) free functions A(p), V(ip) and W{ip) can obviously pass all present tests. 
For instance, GR plus a cosmological constant simply corresponds to A(ip) = 1, V(p) = A/2 
and W{p) = 0. But the presence of a non-constant coupling W{p) can change the physics at 
small scales, notably in the very early universe (Big-Bang) and for later clustering properties. 

The fact that W (ip) induces effects at small scales can be understood by a simple dimensional 
argument. Since this function multiplies the square of the curvature in action (P), it induces 
corrections proportional to 1/r 7 (and higher orders) to the Newtonian potential in 1/r, and 
thereby generically dominates at small scales. However, we saw above that this quick reasoning 
can be erroneous in some perturbative but highly nonlinear situations. Indeed, if Wq takes very 
large and negative values, the cosine involved in Eq. (EJ) shows that p is always of the order 
of —Wq/Wq, even for small distances r. One can then prove that the (very easily satisfied) 
condition | | <C r 2 suffices for all scalar-field effects to be negligible in the solar system, 

even if \Wq\ ~ 10 120 . This remark underlines that nonlinear effects can drastically change 
the intuitive behavior, but let us recall that our cosmological reconstruction above predicted a 
large and positive value for Wq . In that case, we did find that the scalar-Gauss-Bonnet coupling 
induces large effects at small scales, and even exponentially larger than the linear results ©-©. 

5 Conclusions 

Scalar-tensor theories of gravity are the best motivated alternatives to general relativity. Three 
classes of experimental data give qualitatively different constraints on them. Solar-system tests 
strongly constrain the first derivative of the matter-scalar coupling function A(p) (i.e., the 
linear matter-scalar coupling strength). Binary-pulsar data forbid large and negative values 
of its second derivative (quadratic matter-scalar-scalar coupling). The knowledge of the two 
cosmological functions Dl(z) and S m (z) suffices to reconstruct the full shape of both A(<p) and 
the potential V(p) on a finite interval of tp. The knowledge of the luminosity distance Dl(z ) 
alone over a wide redshift interval strongly constrains the theories if one takes into account 
solar-system (and binary-pulsar) data, the positivity of the graviton and scalar energies, and 
the stability and naturalness of the models. Future data, provided by experiments like the SNAP 
satellite, will notably allow us to discriminate between GR plus a cosmological constant and a 
potential-free scalar-tensor theory. The possible coupling W(<p) of the scalar field to the Gauss- 
Bonnet topological invariant can be constrained only if one takes into account cosmological and 
solar-system data together. The predictions of the model at small distances can depend on 
highly nonlinear corrections. Of course, a model including all three functions A{p), V(tp) and 
W{p) is experimentally allowed, since GR plus a cosmological constant is a particular case. The 
presence of a scalar-Gauss-Bonnet coupling W{<p) will generically change the behavior of the 
theory at small scales (clustering, Big Bang). 
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